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Abstract 

We prove a formula for the n-th derivative of the period function T in a 
period annulus of a planar differential system. For n = 1, we obtain Freire, 
GasuU and Guillamon formula for the period's first derivative [17| . We 
apply such a result to hamiltonian systems with separable variables and 
other systems. We give some sufficient conditions for the period function 
of conservative second order O.D.E.'s to be convex. 

Keywords: Period annulus, period function, normalizer, lineariza- 
tion, hamiltonian system, separable variables. 



1 Introduction 

Let be an open connected subset of the real plane. Let us consider a differential 
system 

z'^V{z), z={x,y)en, (1) 

V{z) = {Vi{z),V2{z)) e C^{n,]R^). We denote by </)y(t,z) the local flow 
defined by ([T]). A topological annulus ^ C is said to be a period annulus of 
([T]) if it is the set-theoretical union of concentric non-trivial cycles of Q . If the 
inner component of ^'s boundary is a single point O, then O is said to be a 
center, and the largest connected punctured neighbourhood Nq of O covered 
with non-trivial cycles is said to be its central region. If A is a period annulus, 
we can define on A the period function T by assigning to each point z € A the 
minimum positive period T{z) of the cycle 7 passing through z. We say that 
the period function T is increasing if outer cycles have larger periods. Let 7^(s) 
be a curve of class meeting transversally the cycle 7 at the point s — sq. 
We say that 7 is a critical cycle if [j^T{'-f^{s))]^_^^ — 0. It is possible to prove 
that such a definition does not depend on the particular transversal curve 7^ 
chosen. 

The existence and number of critical cycles affects the number of solutions 
to some boundary value problems. In fact, given a positive r S IR, the number 
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of T-periodic cycles contained in a period annulus A is bounded above by n + 1, 
where n the number of critical cycles contained in A. Similarly, the existence 
of critical orbits is related to the study of some Neumann problems for systems 
equivalent to second order differential equations, as well as to the study of mixed 
problems. The absence of critical orbits is itself an important element in the 
treatment of boundary value problems, bifurcation or perturbation problems 
[3 [32], delay differential equations [TT], thermodynamics |57j, linearizability 
[26]. 

The simplest case, that of monotone period functions, was dealt with in 
several papers. For several references we address to the bibliographies of [T71 
[331 [351 13S] • A very special sub-case is that of isochronous systems, i. e. systems 
having period annuli with constant period function, for which we refer to the 
survey [J , and to the bibliographies of some recent papers as [S] [23] . Systems 
with one or more critical orbits were studied in[2l[^l51[TT 1 [T8 l [T9 l [2T l [27 l [30 l 

In some papers upper bounds to the local or global number of critical or- 
bits are studied [T] |H1 [S] • Finding an upper bound to the number of critical 
cycles is similar, to some extent, to the problem of finding an upper bound 
to the number of limit cycles of a planar differential system. In fact, similar 
techniques have been developped in the treatment of such problems, mainly in 
a bifurcation perspective. In particular, the role played by the displacement 
function's derivatives in limit cycles' bifurcation is similar to that of the pe- 
riod function's derivatives for critical cycles bifurcation [T] [H [12 [T3 [THl US] . 
Dealing with bifurcation from a critical point O, the key piece of information is 
the order of the first non-vanishing derivative at O. On the other hand, global 
estimates to the number of critical orbits are usually based on some property of 
T's derivatives in all of a period annulus. This task has been faced in different 
ways in H [3 [i [H] [H [Ml [ST] for second order ODE's, [S [H [SD] for other 
types of hamiltonian systems, [HI [19] for complex differential equations, [20] 
for polynomial systems. 

In this paper's view, a key result was given in pJTj . Let us denote by [V, W] = 
dyW — dwV the Lie bracket of V and W. A vector field W, transversal to V, is 
said to be a non-trivial normalizer of on a set A C 51 if there exists a function 
IX defined on A such that \V, W] = iiV on A. We call /i its N-cofactor. Let 
(j)w{s, z) be the local flow defined by the solutions of 

z' = W{z). (2) 

In [IT], it was proved that 

dwTiz) = ^ T(0M/(s, z)) = ^ /i(0v'(t, z))dt. (3) 

Hence, if a normalizer is known, such an approach allows to get some information 
about the first derivative of the period function, in particular when the function 
H does not change sign, avoiding the need to evaluate the above integral. In 
order to find a normalizer, it is sufficient to know a first integral, as shown in 
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[50] . In some special cases, more convenient normalizers can be found, as for 
Hamiltonian systems with separable variables [T7]- 

In this paper we give a formula for the n-th derivative of the period function, 
based on an approach similar to that one introduced in |T7j. Rather than 
focusing on a 7V-cofactor, we base our approach on a the existence of a suitable 
commutator, looking for a function m(z) such that [mV, W] = 0. We call m a 
C-factor. C-factors and N-cofactors are related by the equality dwrn = rnfj., so 
that they can be obtained from each other, at least in principle. Let us denote 
by dy^^T the n-th derivative of T{(j)]^(s, z)) with respect to s. In our main 
result, we prove that 

dw = ^ Ticbwis, z)) ^ m{Mt,z)) 

In the case of the first derivative, our formula reduces to Q. We also give a 
recursive formula that, starting from a N-cofactor ^, allows to avoid C-factors. 
Setting 



we prove that 



d\;;^T{z)^ I ^ir,{cj>v{t,z))dt. 







In particular, T is VF-convex if fj.2 — fJ-^ + dwf^ > 0. We also provide a wide 
class of systems having explicit C-factors, including hamiltonian systems with 
separable variables. In such a case, we extend the formula given in |17| for 
hamiltonian systems with separable variables to higher order derivatives. In 
particular, we prove that a hamiltonian system with separable variables 

x' = F\y), y' = ~G'{x) 

has at most one critical orbit if 



^J■s2 = 4 



1 + 2- 



GG" FF" 



G'2 

3G2G"2 - BGG'^G" - G^G'G'" 3F^F"^ - 3FF''^F" - F'^F'F'" 



G'4 F'^ 

has constant sign. As a consequence, we prove that for systems equivalent to 
second order differential equations, the inequality 

G"' - SGG'^G" -I- \2G'^G"^ - 4G2g'G"' > 

implies T's convexity. 
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2 Some properties of normalizers and period func- 
tions 



We assume V and W to be transversal on a period annulus A, and W to be 
a non-trivial normalizer oi V, i. e. V A W 0, and [V, W] — fiV on il. 
\V, W] ■ V 

One has fi = — - — rr — e C°°{A, JR). Let us choose arbitrarily a point z* in A. 
\V\ 

The orbit (fiwis, z*) meets all the cycles of A. Since is a normalizer, the map 
(l>w{s, ■) takes l^-cycles into 1^-cycles, hence there is a one-to-one correspondence 
between V^-cycles and the values of s in some real interval Iw- We parametrize 
such cycles by means of the parameter s. The transversality of V and W implies 
that T G C°°(^,IR). Since, by definition, T is constant on cycles, hence a first 
integral of ([T]), we may write T{s) to denote the period of the unique V^-cycle 
corresponding to the parameter s. Different normalizer can produce the same 
parametrization. In fact, since, by the transversality of V and W every vector 
field definied in A can be expressed as a linear combination Z = aV + f3W, 
/3 > in order to preserve transversality, one has 

[V, Z] = [V, aV + PW] = {|3^JL + dva)V + (9^/3) W. 

Hence aV -\- (3W normalizes V if and only if (9y/3 = 0, i. e. /3 is a first integral 
of ([T]). Moreover, the new N-cofactor is -I- dyce. In a sense, we can split 
the action of aV + /3W on A as the combination of a rotation along ^-orbits, 
determined by the term aV, plus a motion along the VF-orbits determined by 
the term /3W. Applying the main theorem in |17j that gives the first derivative 
of T w. resp. to s, 

dzT = [ {(3fi + dva){(j)v{t, z)) dt = 
Jo 

= I3[ ^li(j)vit,z)) dt = /3 dwT, 
Jo 

we see that the contribution of the term dyct is zero, while the presence of the 
y-first integral /3 only appears as a factor that can be taken out of the integral 
in (|3|, since it is constant with respect to t on (j)y{t,z). In conclusion, every 
normalizer gives the same expression for T'{s), up to a multiplicative factor, 
which is a first integral of ([T]). 

In some cases a particular parametrization for the ^-cycles is preferred. This 
is the case of hamiltonian systems 

x' - Hy, y' - (5) 

where the "natural" parameter is provided by the Hamiltonian function H[z). 
A normalizer producing H{z) as a cycle's parameter is 
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since for such a system one has H = 1. In \WP it was proved that the related 
N-cofactor is 



{Hyy - H^x)H^ - AH^yH^Hy + [H^r, ~ H yy) 

m , (7) 

so that 



T'{H)^ f ' \Hi^v{t,z))dt, (8) 
Jo 

performing the integration afong a T^-cycle (j)v{t,z). For some hamihonian 
systems it is possible to find other normalizers such that H = 1. For instance, 
if H{x,y) — F{y) + G{x), one can consider the system 

,^ G ,^ F 

^ HG'' ^ HF'' 

a special case of the normalizer U* in [T7], remark 7, obtained choosing k{x,y) = 

= —— -. On the other hand, such a normalizer does not exist at points z such 

that G{z) ^ 0, G'{z) = 0. 

If is a non-trivial normalizer of V ^ then it is a non-trivial normalizer of 
pV , for every non-vanishing function p e C°°(ri,IR^). In fact, assuming p > 0, 
one has 

[pV, W] = (pM - dwp)V - (/i - dw (In p)) {pV) = JI ipV). 

This finds an application when a first integral H{x, y) of a system ([I]) is known, 
since the system (|6]) is a normalizer of every system having H(x,y) as a first 
integral. We say that a non-vanishing function p € C°°(Q,TR^) is an inverse 
integrating factor of the system ([T]) if {—V2{z),Vi{z)) — p(z)VH(z), for some 
H{x,y). If this occurs, then the N-cofactor corresponding to the normalizer (|6| 
is 

Ph^ i^H - 9vy(lnp). 

Next lemma shows that differentiating with respect to s produces new first 
integrals of Q. Let us set 

Lemma 1 For every n is a first integral of ([l]) . 

Proof, it is sufficient to observe that both s and T[s) are first integrals of ([T]), 
so that every derivative of the period function w. resp. to s depends only on s, 
i.e. it is constant on the cycles of ([T]). ^ 

In the following we shall denote by 



T^-)[z) = d^-'>T^—T{cj^w{s,z)) 



the n-th derivative of T with respect to s. Dealing with r("), for n > 1, leads 
to change one' approach in relationship to the choice of normalizers. In fact. 
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when studying the sign of T' it makes no difference to use any normahzer, but 
for computational complexity. On the other hand, different normalizers might 
produce higher order derivatives with different signs. Let us denote by a the 
parameter induced by Z, so that (/)n/(s, z) — (f>zicr{s), z), with cr'(s) > 0. One 
has 

^n^wis, z)) = <j\s)^T{c^z{<j{s), z)), 
as as 

while 

f^T{cf>w{s, z)) = a"{s)^T{M<^{s), z)) + {a' {s)f ^T{M<^{s),z)). 
ds'^ as as'' 

This shows that convexity is not normalizer-invariant, so that we shall say that a 

period function T is W -convex, rather than just convex. This opens the problem 

to find the most convenient normalizer W in order to study the sign of d^'^T. 

In general, the unique d'l^h to have normalizer-invariant sign is dwT. If 
a cycle is critical with respect to a normalizer, then it is critical with respect 
to any other normalizer, and a period annulus A can be split into sub-annuli 
where T is monotonic with respect to any normalizer. 

In this paper's applications, we mainly consider the convexity problem and 
its consequence on critical cycle's uniqueness, in particular around a center O. 
This leads to different situations for degenerate centers and for non-degenerate 
ones. In order to illustrate the behaviour of the period function in a neighbour- 
hood of a center, let us restrict to analytic systems. If O is a non-degenerate 
center, then T has an analytic extension at the origin |35) . Since analytic func- 
tions cannot have infinitely many zeroes accumulating at a point, every line 
segment passing through the origin has an open subsegment E such that all 
cycles meeting E have only two points on E. Possibly rotating the axes, we may 
assume E to be contained in the x-axis, i.e. E — {(a;,0) : —e < x < e,e > 0}. 
Then one can define an involution t : E — > H, i. e. a function such that such 
that i'{i{x)) = X, satisfying T{x) = T{l{x)) 24]. For si-reversible centers, such 
a property reduces to T{x) = T{~x), i. e. T is even. As a consequence, its 
Taylor expansion on the a;-axis has the form 

T{x, 0) = a2kx''' + a2k+ix'^^^ + ... 

Hence, if a2k > 0, it is both increasing and convex (decreasing and concave) in 
a neighbourhood of O. Hence, if we want to prove the uniqueness of a critical 
orbit, we can only try to prove convexity/concavity out of a neighbourhood of 
O. If the center is degenerate, then the period function is unbounded at the 
origin, so that T could be both decreasing in a neighbourhood of O and convex. 

3 Results 

We say that a map linearizes a differential system if it takes such a system into 
a linear one. Next lemma was contained in the unpublished preprint [28] as 
theorem 3. 
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Lemma 2 Let A be a period annulus of with \V, W] =0 on A. Then there 
exists a map A e C°°(A,]R2) that linearizes both and 

Proof. By hypothesis, V and W commute, i. e. if (f>Y {t , {(pw {s , z))) and 
4>w{s, {ipvit, z))) exist for ah s and i in a rectangle Is'xit, , then 0y (i, {(pwis, z))) 
= 4'w{s, z))). By the main theorem in 29\, A is an isochronous annulus, 

i. e. every y-cycle in A has the same period T > 0. Possibly multiplying 
V hy we may assume ([T]) to have period 2n. Following 29 , we choose 
arbitrarily a point z* in A and define on A the functions t(z), s{z) such that 
^y(t(z), {(f>w{s{z),z*))) = (l)w{s{z), {(l)v{t{z), z*))) = z. The regularity of t{z), 
s{z) can be proved by the implicit function theorem, as in [29^ or [TJ], section 
4. By construction, one has 

dvs{z) =0 / dws{z) = 1 



dvt{z) =1, \ dwt{z) = 0. 

Let us define A : A i-> IR^ as follows: 

{u,v)=A{z) = (e^W sint(z),e''(^) cosi(z)). 

In order to show that A is injective, first consider that there is a one-to-one 
correspondence between V-cycles and values of s{z). Moreover, there is a one- 
to-one correspondence between points of a y-cycle and the values of t{z), for 
t{z) G [0, 27r). The map A transforms injectively the l^-cycle corresponding to 
s{z) into the circle (e*'^-' sint(z), e*'-^-' cosi(z)). 

Using Q it is immediate to show that A transforms ([T]) and ^ respectively 
into the following systems 



u = V 
v' = —u 



(10) 



In next theorem we give a formula for the n-th derivative of T w. resp. to the 
parametrization induced by a transversal vector field W such that [rnV, W] = 0, 
for some non-vanishing multiplier m. We may assume m to be positive, since 
the period function of V and that one of —V coincide. 

Theorem 1 Let A be a period annulus of ^ and m G C°°{A, IR), m ^ 0, such 
that [mV, W]^0 on A. Then 



" Jo m(0r(t.«)) 

Proof, by hypothesis, the systems 

z' = m{z)V{z) z' = W{z) 



dt, (11) 
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commute. By lemma [2j there exists a transformation A G C°°(A,]R^) which 
takes such systems resp. into the systems (10). As a consequence, the system 
([TJ is transformed into the system 

( u' ^ - 

^™*(«^ , (12) 



where m*{u,v) — m{A~^{u,v)). Since A preserves the time, the system (12| 
has a period annulus with cycles having the same periods as their anti-images 
in A. Denoting by 0* the argument function in the plane (m, u), one has 

r= / ldt= ^ = - m*dd*. (13) 
Jo Jo 0* Jo 

The system ^2]) is transformed into the second system in (10), which is a nor 



malizer of (12). We denote its vector field by W* . We can find the derivative 



of T with respect to the parameter s by differentiating the integral in ( 13 ) with 
respect to s: 

dw'T = = / m*de* = - / ^m*d6* ^ / — ^m* dt. (14) 

ds ds Jq Jo ds Jq to* ds 

As for higher order derivatives, a similar conclusion holds, since the integration 
extremes do not depend on s: 



Ou/.i = -, — 1 — — / m dO ~ — I -r — TO du — -r — TO dt. 

^ ds"" ds'"- Jq Jo as" Jo to* 9s" 

(15) 

Since A preserves the time of all the systems considered, applying the inverse 



transformation A gives the formula (11). Jfk 



In next lemma we show that the relationship [mV, W] = 0, for a non- 
vanishing TO, is equivalent to W being a normalizer of V, and find the rela- 
tionship between to and the N-cofactor /i. We state it for a period annulus, 
even if it holds in other subsets of il. 

Lemmas Let V,W e C°°{A,M^). Then 

i) there exists to G C°°{A, IR), to(z) > 0, such that [mV, W] — 0, if and only 
if W is a normalizer of V , with cofactor fi = 5vi^(lnTO); 

a) m,m e C°°{A,]R), m{z),m{z) > in A, satisfy [mV,W] = = [toV, W], 
if and only if there exists J € C°°{A,1R), first integral of such that 
fn — Jm. Moreover, 

dw"^ ^ dw'^ _ (16) 
to to 
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Proof, i) If [mV, W] = 0, then 

= [mV, W] = m[V, W] - {dwm)V, 

hence 

[V,W] = (^^^ V = (^dw{lnm)y. 

Conversely, let us assume there exists /i such that [V, W] — fiV. Let us choose 
arbitrarily a point z (£ A. Let t{z), s{z) be the functions defined in in lemma [2j 
so that z = (f)v{t{z), (j)wis{z),z)) = 4>wis{z),(f)vit{z),z)). Let us define m{z), 
as follows: 

One has m e C°°(yl,IR) and 

dwm{z) = ^,(^)e/;'^VWw(<.,0v(t(.),^))d<. ^ ^(z)m(z), 
since dws{z) — 1. Then 

[mV, W] = m[V, W] - {dwm)V = mfiV - fxmV = 0, 
ii) Let J be a first integral of Then 

[JmV,W] ^ Jm[V,W]-{dw{Jm))V = j(mfi-dwm^V ^ J {m^iV - ^imV) = 0. 
Vice- versa, if [mV, W] = = \jriV. W] , then one has 

(m/i — dwrnjV = [niV, T^] = = [mV, W] ~ (rfifi — dwrn)V. 
From n = = one has 

" m m 

m dwrfi — m dwm — 0, 

that gives 



TO 



TTl 

hence — is a first integral of ( 2 1 . 

TO 1_| 

As for (16), if TO = Jto, then d'^^ m — (Jm) — J d'^^ to, since J is a 
first intcgraJof ([2|. Hence the two fractions in pG] ) coincide. Jfk 

In lemma|3]the integration is performed only along the M^-orbits, starting at 
the intersection of the ^-cycle through z with the VF-orbit passing through z. 
Performing the integration along the W-orbit, we have chosen an initial value 
of 1 on the cycle so that m{(j)v{t,z)) = 1 for all i € H. We could 

have chosen other smooth initial values on getting other isochronous 

systems with the same period annulus A. All such functions to can be obtained 
from one another multiplying by a first integral of W. All of them give the same 
ratio that appears in the integral ( [Tl| ). 

As a special case of lemma [3] we obtain the cited formula for the first deriva- 
tive of the period function [17] . 
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Corollary 1 Let A be a period annulus of with \V, W\ = on A. Then 



dwT{z) — J fi{<j)vit, z))dt. 
Proof. It is sufficient to observe that 



n ^ 9 r * ^ fj,(cl>iy(ij,<pv(t(z),z))da / % f ' ^ ii(4>w{c,'l>v(t(z),z))d(j 

OS OS 

hence 



^— = [i{z). 

m[z) 



In some cases it is computationally more convenient to look for a fi such 
that [V, W] — nV, rather than for a m such that [mV, W] = 0. In next theorem 
we give a direct way to compute high-order derivatives of T working only on /i. 
This can be done by means of a recursive formula. 



Theorem 2 Let A be a period annulus of with [V, W] = fiV on A. Th. 



en 



d^w^T{z)= r t,^{<j,v{t.z))dt. 
Jo 

where /z„ is recursively defined by 

Ml = /"n = A*n-lM + dwlJ-n-l- (17) 

Proof. Let us define m as in the proof of lemma [3j Then, let us set 

m„ = . 

m 

By corollary [TJ one has rrii = fi. For n > 1, one has 



/gin-l 


TO^ 


\ m 








TO 


TO 



= TO„ - TO„_lTOl, 

;/t //t TO 

hence 

TO„ = TO„_iTOi + dwrrin-i. 



This shows that the functions to„ satisfy the recurrence equations (17), hence 
TOrj = Hn- Then the statement comes from theorem [l]. ^ 
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In the next formulae we replace the notation dwl^ with the simpler one /z'. 
Similarly for higher order derivatives with respect to W. We write the forms of 
some low order /x„'s: 

^4 = / + e^iV + 4/iM" + + /i'". 

= ^5 + 10^3^y ^ 10/xV" + 15AiM'^ + 5/i/i'" + 10/iV" + ^i^'^^ 

We say that a fmiction /i G C°°(A,IR) satisifies the condition (B) if /i(z) 
does not change sign in A, and every cycle contained in A contains a point z 
such that h{z) ^ 0. 

Corollary 2 Let A be a period annulus of with [mV, W] = 0, with m > 0, 
and [V, W] = /iV^ on A. // one of the following holds, 

i) dy^'^m satisfies the condition {B), 
ii) fXn satisfies the condition (B), 

then A contains at most n — 1 critical cycles. 

Proof, i) T is a first integral, hence every critical cycle of ([l]) corresponds to 
a critical point of T{(j)w{s,z)). By ([ll|, if d\^'^m > 0, then d^^'^T > 0. The 
presence on every cycle of a point where d^^T ^ implies that on every cycle 
dl^'^T > 0. Similarly if d[^''m < 0. 

ii) If fin satisfies the condition (B), then also d^m satisfies the condition 
(B). 4k 

In particular, the corollary [2] holds for n = 2, in presence of a convex or 
concave (w. resp. to the parametrization induced by VF) period function. The 
convexity may be proved under a condition on dwfi, rather than on d^J m. 

Corollary 3 Let A he a period annulus of with [V, W] = fiV on A. If 

dwf' > on A then T is W -convex. If dwfJ- satisfies the condition (B), then T 
is strictly W -convex on A. 

Proof. By theorem [2] one has 

ti2 — f-^ + dwf- > 0. 

Then 

a{^^r(z)= / n2{<f>vit,z))dt>o. 

Jo 

If every cycle contains a point where fi2 > 0, the above integral is positive. Jit 

Remark 1 It is remarkable the asymmetry of such a situation, where convexity 
can be proved by only considering the sign of dwf^, while concavity cannot. This 
agrees with the fact that T can be upper unbounded, hut it cannot be lower 
unbounded. 
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4 Applications 



4.1 Reparametrized isochronous centers 

A center is isochronous if and only if it has a non-trivial commutator |29j . 
Given an isochronous center, finding a commutator is not always easy. There 
are systems which are known to have isochronous centers, for which no com- 
mutators are known, as for reversible Lienard systems [2 HUl H] . A collection 
of isochronous systems with their commutators have been provided in [1] , other 
isochronous systems may be found in [31 HHj and their bibliographies. 
In this section we consider centers of the form 



(18) 



where V satisfies [V, W] = on a period annulus A, W transversal to V, > 0. 
By lemma 3 an obvious choice to study the period function of is m — -. 
Hence we have the following corollary. 



Corollary 4 Let A be a period annulus of CT, and W be such that [V, W] = 0. 
1 J I 

// - satisfies the condition (B), then the system (18) has at most n critical 

cycles. 



Proof. It is an immediate consequence of corollary 



choosing m 



It is not difficult to find examples of systems ( 18 ) with exactly n critical 
cycles. It is sufficient to consider a linear center 



x'^y^{x,y), y' ^ -x£,{x,y), 



(19) 



with £^{z) — uj{\z\), ^jlyj having exactly n critical points. In fact, in this case 
one has T{z) = = 

Replacing again the notation dwlJ- with /i' one has 



dwm _ „ 

Ml = M = — =iOw 
m 



M2 = M 



/i^ + + /i" 



2e - c'^ 
e ■ 

-6^'^ + eee'c" 



/X4 = /-l-6/iV+4/x/x"+3^'2+/ 



24C^ - secede + s^^cc" + 6^^^"^ - C^C^^^ 



/i5 = /i''' + 10/i'^/i' + 10/i'^/i" + I5fifi'^ 



5fj.fl" 



IQfi fi' 
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Let us assume i{x,y) to be analytic, ^{x,y) — J2'^=o^n{x,y), where ^„ = 
X]J=o '^n,jx"'~^y^ is an n-degree homogeneous polynomial. If ^ > in a neigh- 
bourhood of O, then the origin is a center of (19). One has V{x,y) — (y, — x), 
W{x,y) = {x,y), so that, by Euler's formula, 

(oo \ / oo \ oo 

^^n{x,y)\ +y i^C7i{x,y)\ =^n$„(x,y). 
n=0 ) X \n=0 / y n=0 

Similarly 

oo 



Then, working as in corollary |3j if 

oo 



n=0 



then T is Ty-convex, since /i2 = /i^ + /i' = — > 0. Moreover, if ^ 



e 

satisfies the condition B, then T is strictly ly-convcx. This is the case of the 
system 

a;' = -£,n{x,y)), y' = -a; (1 - ^„(x, y)), 

if £,n{x,y) > 0, ^nix,y) homogeneous of degree n. On the other hand, since 
— ^' = n^nix,y) > 0, T is increasing, hence there are no critical orbits. In 
this case the period annulus is contained in the oval defined by the inequality 
£.n{x,y) < 1. 

In next corollary we consider a class of systems admitting critical orbits. 

Corollary 5 Let ^i{x, y) homogeneous of degree I for I = k,n. Assume k < n, 
(3 - 2V2)n < /c < (3 + 2V2)n, ^k{x,y) > for {x,y) ^ (0,0). Then O is a 
center of 

x' = y{£,k{x,y) -^n{x,y)), y' = -x{£,k{x,y) -^n{x,y)), 

with period function T W -convex on Nq- If, additionally, ^n{x,y) < 0, then 
there exists exactly one critical orbit in Nq ■ 

Proof. The origin is a center, since the condition > implies that in a 
neighbourhood of O the orbits coincide with those of the linear center. The 
numerator of 

/i2 = yU^ + fl' = 



2 e - 
e 



IS 
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Considering it as a quadratic form in tlie indeterminates f^, a sufficient 
condition for /i2 not to cliange sign is 

A = (fc2 + - iknf - ik'^n^ = {n^ - 6n/c + k^){n - kf < 0. 

If (3 - 2V2)n <k<{3 + 2^/2)n, then - Qnk + fc^ < 0, hence A < 0. 

If, S,nix,y) < 0, then the central region Nq is contained in the oval having 
equation £.kix,y) — £,nix,y) ~ 0. Such an oval consists of critical points. The 
system is analytical, hence the boundary dNo cannot be a limit cycle, so that 
at least one of such critical points lies on the external boundary of Nq- As a 
consequence, T goes to oo approaching the oval. Moreover, T goes to oo as well 
approaching O, since the center is degenerate. Hence T has a minimum in A, 
reached on a critical orbit (j)v{t,z), which is the unique critical orbit in A by 
T's W-convexity. ^ 



4.2 Jacobian maps and Hamiltonian systems with separa- 
ble variables 

Let ^{z) = (P(z),Q(z)) € C°°{n,]R^), with jacobian matrix J^,{z). li 5{z) = 
det J^(z) 7^ 0, we say that it is a jacobian map. If ^I' is a jacobian map, 
possibly exchanging P and Q, we may assume its jacobian determinant S{z) = 
det J*(z) to be positive on SI. This ensures local, but not global invertibility. 
Let us consider the function H{z) defined by 2H{z) ~ The Hamiltonian 

system having H as Hamiltonian is 



Hy^PPy + QQy, v' = - H , = - P P, - Q Q , . 



We consider also the system obtained dividing ([T]) by 5{z), 



PPy + QQy 



(20) 



(21) 



and the system 



y = 



-PQ, + QP, 



(22) 



From now on, we denote by the vector field of Q, V5 the vector field of (21 ), 
W5 the vector field of (22 1. We denote by (t)y^{t,z) a solution to ([T]), and by 
4>Ws{^tz) a solution to (22 1. 

In [3S] critical points of (21) were proved to be isochronous centers, under 



the assumption (5 > 0. In fact, the systems (21 1 and (22 1 commute with each 
other, that allows us to apply the theorem ijto find a formula for the n-th 
derivative of some hamiltonian period functions. We denote by d^^^T the n-th 
derivative of T with respect to the function H. 
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Theorem 3 Let ^{z) = (P(z),Q(z)) G C°°{n,M'^). If A is a period annulus 
of p^, then 



dt = 



T(z) 



(23) 
(24) 



P^iin{4'V^{t,z))dt, 

where fx^n is defined recursively as in theorem^ with ji^ — — In (5. Moreover. 



Jo 



lJ'n{4)V^{t^z))dt, 



(25) 



where /i„ is defined recursively as in theorem 



2: with fi 



In (5 
2H' 



Proof. The transformation 5* takes locally the systems (21) and (22 1 into the 



from (21 



systems ( 10 ), which commute with each other. Since the system i20 ) is obtained 

. . . 1 

multiplying by S, setting m = - one has [mVg,, W^] = 0. Then the 





equality ( 23 ) comes from theorem [Tj By lemma |3j one has 



A** 



5 ^^ 



(52 



-(9vK5(ln^). 



The formula ( 24 1 conies from theorem [2] 

The function H vanishes only at critical points of V*, which coincide with 

those of Yi, and W5, hence the vector field — — - is a non-trivial normalizer of 

2i/_ 

ln(5 



20 ) , with N-cofactor /i = . The formula (|25|) comes as well from theorem 

w 

In theorem's |3] proof one does not need the invertibility of ^ on all of the 

period annulus, since one only needs a C-factor, provided by -. 

1-1 . . / " . 

The applicability of theorem ^\ depends on the possibility to write a hamil- 

tonian function ii{z) in the form 2H{z) = \'^{z)\'^, with det J*(z) ^ 0. Such a 
question was addressed in [33l |26j . 

Since the above theorem can be stated in terms of normalizers, we can con- 
sider also non-hamiltonian systems, as observed in section 2. 

Corollary 6 Let p > be an inverse integrating factor of the system with 
first integral H(x,y) satisfying 2H{z) = |^(z)|2j for some jacobian map '5. // 
A is a period annulus of (20), then 



rT(z) 

d\;^lTiz)^ / fini^vAt,z))dt, 
Jo 



(26) 
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where /i„ is defined recursively as in theorem^ with /i = —\n{5p). Moreover, 

fT{z) 

d^H^T{z)= / Ti,XMt.^))dt, (27) 
Jo 



where /i„ is defined recursively as in theorem 



2: with ji 



ln((5p) 
2H 



Proof. By hypothesis, one has {—V2(z),Vi{z)) = p{z)WH{z), with 2H(z) 
|^(z)P, hence V{z) — p{z)V<s,{z). Working as in section 2, one has 

[V, Ws] = [pV^, Ws] - - dw,p)V = 



dwAl^S)~dws{^rip) UpV) 



[-dwsi HSp))){pV^) = ( - dws ( H5p)))v. 



that gives the formula ( |26[ ). The formula (27 1 can be proved as formula (25 1 in 
theorem [3] ^ 

In this paper we only consider the applicability of theorem [3] to some hamil- 
tonian systems with separable variables. Let Ip, Iq be intervals containing 0, 
F e C°°(/f,IR), G e C°°(/g, IR). Assume G{x) and to have isolated 

minima at 0. Then the origin O is a center of the hamiltonian system having 
H{x, y) ~ G{x) + F[y) as hamiltonian function, 



F'iy) 
-G'{x). 



(28) 



Un der quite general conditions, such systems can be considered as special cases 
of (|20|, obtained taking -^(x.y) = {P[x),Q{v)) = {s{x)^2G{x), s{y)^2F{y)), 
where s{t) the sign function, assuming values —1,0, 1 for t < 0, t — 0, t > 0, 
respectively. The jacobian determinant of such a map is S{x, y) — P' {x)Q' [y] = 



F'(v)G'(x) 

s{xy) — 7^^=^^^^^- In this case, the normalizer (|22|) has the form 



2v/F(2/)G(a;)' 

PQy - QPy ^I^^2G 

S ^ P' ^ G'' 



y 



-PQ. + QPx ^ (I 
S Q' 



2F 



(29) 



and differs from that one given in [T7j for the presence of a factor 2. Also the 
corresponding N-cofactor 

Psix,y) = -avK5(ln(5(x,?/)) = -dws[MP'ix)Q'{y)j ^ 



In 



s{x)G'{x) 
v/2GM 



In 



<v)F'{y) \ ^ (G(x)\ (F{y)y 



differs from that one given in |17j for the presence of a factor 2. The commuting 
system (21 1 has the form 

, _ 2s{xy)^G{x)F{y) _ Q{y) , _ 2s{xy)^ G{x)F{y) 



G'ix) 



p'{xy 



F'iy) 



Pixl 

Q'ivY 

(30) 
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Let us set 



1 



2s{xy)^G{x)Fiy) 



6{x,y) G'{x)F'{y) ' 

whenever F{y)G{x) > 0, F'{y)G'{x) ^ 0. We say that a function i : ^ R, 
Ih open interval containing 0, satisfies the hypothesis (Sep) if 

sepi) L £ C°°(/l,IR), i(0) = 0, iL'(i) > in /l, L is not flat at 0; 



The condition sepi) implies L{t) > for t close to 0, t 7^ 0, hence its Taylor 
expansion starts with an even power of t, 



As a consequence, 
so that 



L(i) = C2fet2'= + o(<2fe). 



4fc2c2^,i4fe-2 + o(<4fe-2) ' 



Such a function can be differentiable only if 4fc — 2 < 2fc, i. e. fc < 1, hence 
A; = 1. This agrees with the fact that non-degeneracy is a necessary condition 
for a center's linearizability j35j . 

Theorem 4 Let F and G satisfy {Sep) on Ip, Ig, open intervals containing 0. 
Then 



i) the transformation ^{x,y) — (^s{x)^j2G{x), s(?/)-\/2i^(y) 
systems and (30); 



linearizes both 



ii) O is an isochronous center of the system (30); 



Hi) the period function of the system (28) satisfies 



ms{(l>v^{t,z)) 



dt, 



(31) 



where the integration is performed along the cycle (f>v^it,z) starting at z. 



Proof. One can prove i) by direct computation. The systems (30) and (29) are 
transformed into the systems (10). 

In order to prove ii), it is sufficient to observe that the systems ( 10 1 commute, 
hence also (30) and (29) commute. By the main result in [29], every period 
annulus of the system (30) is isochronous. 

The statement Hi) is an immediate consequence of theorem [sj Jft 
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The theorem |4] does not apply to degenerate centers of hamiltonian systems. 
In order to study such systems' higher order derivatives, one can apply the 
corollary [2] starting with a suitable normalizer. 

A "universal" normalizer, as pointed out in section 2 (see [31]), that allows 
to cover both centers and period annuli, is 



which gives the N-cofactor 

_ iF''-G"){G''~F'') 

This allows to cover both degenerate and non-degenerate hamiltonians, obtain- 
ing T's derivatives without need to explicitly find the corresponding to, by using 
the recursive formulae of theorem [2j Moreover, such derivatives are computed 
with respect to the variable H . On the other hand, even in the simplest non- 
degenerate cases, this leads to more complex computations than using TOc, and 

the related normalizer fis- As an example. If H{x, y) — , one has 



a;2(3 + 2x2) 



Qx'^ix'^ + Ax-^ + Ax^ - y2) 



The same N-cofactor can be obtained by using the system 

G F 



HG'' ^ HF'' 

which is defined on all of the central region Nq , but is not defined in a neigh- 
bourhood of every point {x, y) where G'{x) = and G{x) 7^ 0, or F'{y) = and 
F{y) ^ 0. This is the case of the hamiltonian system 

x' = y, y' = —X + 2x^ — x^. 

One has 

_ x{3x^ -9x + 8) 



12(1 -x)3 ' 

which diverges at 1. The hamiltonian system has a center at O, with central 

f 2x^ x^ y'^ 11 

region defined by Nn = < (x , y) : < 1 1 < — >. The central 

6 .y u Y 2 3 4 2 12/ 

region's external boundary consist of a critical point at (1, 0) and a homoclinic 

orbit having (1,0) as a- and w-limit set. Every orbit external to Nq is a cycle 

enclosing Nq- One cannot study T's derivatives on such cycles by means of 

Us, since every external cycle has a point (actually, two points) where /i^ is not 
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defined. Obviously, also the jacobian map approach can no longer be applied, 
since 6 vanishes where G' vanishes. On the other hand, we can compute for 
the normalizer 



(x - 2a;2 + 2:3)2 + y2 ' V [x - 2x'^ + x^y + y'^ ' 

_ x(-4 + 3x){y + X - + x^){-y + x- 2x'^ + x^) 
{x — 2a;2 -|- x^)^ + y2 



y 



One has 



The study if T's monotonicity requires an additional study of jiH, in order 
to determine its sign. Higher order derivatives require the study of much more 
complex functions. 



When dealing with the period function of an analytic system ( 28 ) in a central 
region Nq, the second, and more convenient approach, consists in applying the 
theorem [2] to fis- Assume 

G{x) = 0(^2'=), F{y) - 0{y''^). 



Then the system (29) is of class C°° on Nq, as well as /i^. Applying the theorem 
[2] does not require to produce the corresponding C-factor, even if in some simple 
cases it can be found, as for the systems 



X ^y 



The normalizer ( 29 ) is 



X 



„2fe-l 



y 



fee IN, A: > 0. 



y 
fc' 



4(fc - 1)2 

— . There exist infinitely 



2(fc-l) 2 

with N-cofactor /i^ — , and fj,s2 = Ms 

k 

many C-factors, all producing fj,s as N-cofactor. Among them, one has all the 

functions m(x,y) — — -, with v homogeneous of degree 2fc — 2. 

v[x,y) 

In order to perform next computations with the N-cofactor /is, we find 

convenient to write next formulae in terms of the functions r(a;) = \~\ ^ 

G'(x) 



2F(y) 

^(y) = — using the normalizer (|29l 
F'{y) 



r{x), y' = $(2/), 



as well as its N-cofactor 



^is{x,y)=T'{x) + ^'{y)-2. 
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One can write fig as follows, choosing a and /3 such that a + /3 ^ 1. 

(2 - a)G"2 - 2GG" (2 - - 2FF" 



fJ-s 



(r'(x) - a) + (a>'(y) - /3) 



G"' 



^'2 



Hence, if there exists a couple (a, /3) such that on a period annulus (2 — a)G'^ — 
2GG" > 0, and (2 - l3)F''^ - 2FF" > 0, then T is increasing. For systems 
equivalent to second order conservative equations one does not have such a 
freedom of choice, since 



ny) 



r 

2 



<f'{y) 



f 2G{x) 

[g^ 



G'2 - 2GG" 
G^ ■ 



As a consequence, if G'^ — 2GG" > in a period annulus A, then T is increasing 
in A (see [IT]). 



Let us denote by Wg the normalizer vector field ( 29 1 . One can compute 

Ms2 = A^s + dw.^'S- 

J, „GG" i^F" 
A*^2 = 4^1 + 2— + 

3G2G"2 - 3GG'2G" - G^G'G'" SF^F"^ - SFF'^F'' - F'^F'F'" 



G'4 



Such a formula has been written in such a way to emphasize the presence of a 
single term with mixed variables. 

If O is a center of (28), we denote by TVg, Nq the projections of Nq on the 
X-axis and on the y-axis, respectively. 

Corollary 7 Assume G{x) = h2kX^^ + rg{x) £ C°°{Ig,M), F{y) = C2hy'^^ + 
rf{y) g G°°(/f,]R), e /g n /f, < e in, 62fc,C2/, > 0, r<,(x) = 

o{x^^), rf{y) = o(y2''). Then O is a center, S is of class C°°{No,lR^) , 
Us e G°°(A^o,ni^)- //M52 > ('< O; zn A^g X N^then T is Ws-convex (con- 
cave ) on No ■ If, additionally, on every cycle of Nq there exists a point where 
such inequality holds strictly, then T is strictly Wg-convex (concave). 

Proof. The functions F'{y) and G'{x) vanish only at 0. As for the regularity of 
(29) and fig at x = and y = 0, one has 



G'(x) 



'32k 



2kb2 



52fe -r ' g\ 

which is the product of two functions of class G°° at O. Similarly for F{y). 
Since /i„ is obtained multiplying /I's derivatives by the components of (29), one 
has /i„ S C°°. Then the (strict) IV^-convexity comes from the formula 

rTiz) 



dl^lT{z)^ f fis2{cl>vAt,z)) dt. 
Jo 
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Looking for simpler conditions which imply T's convexity, one can write 

Ms2 = (r' + $' - 2)2 + r"r + $"$ = 
= (r' + $')^ - 4r' - 4$' + 4 + r"r + 

Here appears in two ways as the sum of a squared term plus terms depending 
only on x or y. This allows to study its sign in a simpler way than with genuine 
two- variables functions. We say that a function L{t) satisfies the condition {C\) 
if 

Ca = 2L2i"2 + LV'^L" - L^L'L'" > \L'\ A G H. 
We say that a function L{t) satisfies the condition (C) if 

C = 2iL"2 - L'^L" - LL'L'" > 0. 

Corollary 8 Under the hypotheses of corollary^ if one of the following con- 
dition holds, then T is Ws-convex. 

i) G satisfies (Ca) on Nq, F satisfies (C/j) Nq, respectively, with a + l3 — 3 > 

0; 

a) both G and F satisfy the condition (C) on Nq and Nq, respectively; 
Proof, i) One has 

i^s2i = -4r' - 4$' + 4 + r"r + = 

"2G2G"'2 + GG'^G" - G^G'G'" 2F'^F"^ + FF'^F" - F'^F'F'" 



G'i F'^ 

> 4(a + /3 - 3) > 0. 

hence ^Ji,2 = (r' + ^'f - AT' - 4$' + 4 + T'T + $"$ > 0. 
ii) One has 



-3 



> 



^^s2^^ = r"r+$"$ = 4 

Since 



2G2G"2 - GG'^G" - G^G'G'" 2F'^F"'^ ~~ pF'^F" - F^F'F'"' 



+ - 



Q,4 p,4 



2g2g//2 _ GC'^G" - G'^G'G'" _ G_ 
G^ ~ (74 



2GG"2 - G'^G" - GG'G" 



and G > in TVo, under the the condition (G) one has T'T > in A'^o- The 
same holds for hence ^is2 = (r' + $' - 2)^ + r'T + > 0. ^Ik 



An important special case is that of second order conservative ODE's. In 

s case F{y) — —. In next corollar) 
corollary [8] specialized to such a case. 



2/2 

this case F{y) = — . In next corollary we just write the conditions i) and ii) of 
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Corollary 9 Under the hypotheses of corollary^ if F{y) = one has 
G'^ - m'^GG" + 12G^G"^ - AG^G"'G' 

Moreover, the condition i) of corollary^ holds if and only if 



tis2 



GG'^G" 



G^G'G'" 



2G'^ > 0. 



The condition ii) of corollary^ does not change. 
We can apply the corollary [S] ii) to the system 

x' = 2sinj/cosy, y' = — 2sina:cosa;, 



(33) 



obtained taking G{x) = (sinx)^, F{y) = (sinj/)^. The origin is a center whose 
central region is a punctured open square having vertices at (f,0), (O, f), 
(-f ,0), (0,-f). One has r(a;) = tanx, r"(.T)r(.T) = 2(tan2 a;)(l + tan^ x) > 
in No- Similarly for $(y). 



Figure 1: The central region of system (33 1 




We have an example of system satisfying the hypotheses of corollary|9]taking 



G(x) = e('=°"'^~3'=°"^\ F{y) = y. The related system is 



y 



-3 sin xe 



(cos a;— 3 COS a:) 



(34) 



The origin is a degenerate center, with bounded central region Nq- Its bound- 
ary dNo contains a critical point, hence T diverges as a cycle approaches 

1 1-1-3 cos^ r 

dNo (see Figure 1). One has r{x) = ^ , 3 , T"{x)r{x) 



3* o ^ 
sm X 



3 sin X 



>0, 



$"(z/)$(?/) = in No, hence T is strictly convex. This implies the existence of 
a single critical point, since T goes to 00 both as a cycle approaches O and as 
it approaches dNo- 
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Figure 2: The central region of system (34 1 




The possibiUty to consider separately the terms depending only on x and 
those depending only on y implies that we may swap the components of the 
above systems, obtaining a new system with convex period function, 

x' = 2 sin 2/ cos y, y' - -3 sin^ xe^™^' ""^ . 

In [211 [30], potential functions of the form 
consider here some other rational functions. 



were considered. We 



ax" 



Corollary 10 Let G{x) = 



P{x) 



-, with P{x) = ax^ + hx^ + cx^. // a > 0, 

,6 1 01 0^31,2 



\ + P[x)'' 

6 > 0, c > 0, ahc 7^ 0, 20a''c < 2^1\?ac + 1940^02 + 186'* + 90a« + 219a^6^ 
Sa'^c < 40a^ + 99a^6^, then O is a center, and T has exactly one critical orbit 
in No- 



Proof. One has G'{x) 



P'ix) 

. In a neighbourhood of O, xG'(x) > 0, 

(1 + P(x))^ 

G'{x) = if and only if x = 0, hence O is a center. If G'{x) = only at 
a: = 0, the central region Nq is a strip defined by y^ < 2, since lim G{x) = 1. 

a;— v±oo 

The vector field is bounded on such a strip, hence T goes to +oo as the cycle 
approaches the boundary ONq- The same occurs as the cycle approaches the 
origin, since O is a degenerate critical point. 

If there exists a point x ^ such that G'{x) = 0, then the central region has 
at least a critical point on its boundary, hence T goes to infinity as the cycle 
approaches dNo- Moreover, one has 



1 



Dix) 



36x^^c^ + 



x^%c^ + (66ac^ + lAh^c) 
(48a62 - I2c^ + 400^0) x^ {ila% - 196c) x^+ 



110a6c + 2163) a;i" 
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+ (lOa^ - 66^ - 18ac) x'^ - %ahx^ - 2c? 

with D(x) = (2a + Zbx^ + /^cx^f. Finally, 

_ A{x)^B{x) 

with = (53462 ac + 388a2c2-40a''c + 366'^ + 180aS+438a362)2;8 + (_8a3c_^ 
40a^ + 99a^b^)x'^, B{x) polynomial of degree 32 with positive coefhcients, such 
that -B(O) — 4a^ > 0. Under the condition given in the hypothesis, A is non- 
negative, hence T is strictly convex. ^ 

Let su consider the following functions: 

i) the functions — ; ^, with a, 6, c > 0, hlP' < 56ac: 

jj) the functions —rrr with k positive integer, a,b,c > 0, 2b^k < 

1 + bt-^" + ct°'^ 

+m8a'^ck + 18a'^c + b^. 

Corollary 11 LetG{x) and F{y) be of type j) orjj). Then the period function 
of the corresponding system (2^ is convex. 



Proof. Under the given hypotheses, all the functions of the points j), jj), jjj) 
satisfy the condition (C). In fact, computing the expression ( = 2LL"^ — L'^L" — 
LL'L'" for such functions gives the following. 

, 64a;4(6&c2xi2 + (SGac^ - 5b^c)x^ + ISabcx'^ + ia9) „ , 

[a + bx^ + cx^f >0,rfa,5,c>0, 

56^ < 56ac; 
jj) 

A{x) + B{x) 



2 



D{x) 

with D{x) — x'^{a + bx^ + cx^Y , B{x) polynomial of degree 32fc — 2 with pos- 
itive coefficients, A{x) = {-2ab^k + lOSa^cfc -I- 18a^c -I- ab^)x^^''~'^ + (54fc - 
9)ac^x^^^~'^ + {l&k — ^)bc^x^'^^~'^ . Since A: is a positive integer, one has 54A: — 9 > 
0, 16A:— 4 > 0. The same holds for the coefficient of x^^*''"^, under the hypothesis 
2b^k < +108a^ck + ISa^c + b^. 

By corollary [sj ii), if both G and F satisfy one of the above, one has T's 
convexity. A 

P(x) 

We cannot include in the above list the functions ^— of corollary 

1 + Fix) ^ 

since in that case one has 

n^2^ii2 r<n''in" n'^ni n'" — —60^6 + o(l) 

(1 + ax'^ + bx^ + cx^Y> 
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